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ABSTRACT 

Using well known Lagrangean techniques for uncovering the gauge 
symmetries of a Lagrangean, we derive the transformation laws 
for the phase space variables corresponding to local symmetries 
of the Hamilton equations of motion. These transformation laws 
are shown to coincide with those derived by Hamiltonian meth- 
ods based on the Dirac conjecture. The connection between the 
Lagrangean and Hamiltonian approach is illustrated for first class 
systems involving one primary constraint. 
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It has been conjectured a long time ago by Dirac [1] that local symmetries 
of the Hamilton equations of motion are intimately tied to the existence of 
first class constraints, and that the generators for the infinitessimal symme- 
try transformations are given by the full set of first class constraints. The 
problem of confirming this conjecture has been addressed by several authors 
[2-8]. In this paper we shall present an alternative approach to this problem, 
based on purely Lagrangean techniques. 

On the Lagrangean level there exists a well known method for uncovering 
the local symmetries of a Lagrangean [9,10]. This method leads to gauge 
transformations which are parametrized by a set of arbitrary functions of 
time, whose number equals the number of so called "gauge identities" gener- 
ated by the Lagrangean algorithm. In this letter we shall make use of purely 
Lagrangean methods to derive the transformation laws for the coordinates, 
momenta, and Lagrange multipliers, which correspond to local symmetries 
of the Hamilton equations of motion. In particular we will show that the 
transformation laws extracted from the gauge identities are precisely those 
derived within the Hamiltonian formalism based on the Dirac conjecture. 
The merit of our approach is that it is simple, and that its connection with 
Hamiltonian methods based on the Dirac conjecture is very transparent. The 
basic ingredients going into the proof will be illustrated for the case of first 
class systems exhibiting one primary constraint. 

Consider a system whose dynamics is described by a "singular" Lagrangean 
L(q, q) leading to a set of primary constraints {Qe}. Such constraints, follow- 
ing alone from the definition of the canonical momenta, must necessarily be 
present if the theory is to exhibit a local symmetry. As shown by Dirac [1] 
the corresponding Hamilton equations of motion are given by (i — 1, 2, • • •, n) 

an ea, 

dqi Y d< li 
together with the primary constraints 

n t (q,p) = 0, (3) 

where H(q,p) is the canonical Hamiltonian evaluated on the primary surface 
defined by (3), and where the t/'s are undetermined velocities playing the role 
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of Lagrange multipliers. These equations follow from a variation principle 
5 St = 0, where 



S T = dt 



,i=i e 



(4) 



is the "total" action, and the variations are taken with respect to qi, pi and 
v. 

We now notice that equations (1) and (2) are identical with the Euler- 
Lagrange equations of motion derived from the first order LagrangeanP] 



2n 



HQ, Q) = £ a a (Q)Q a - H T (Q) , (5) 



a=l 



with Ht, the "total Hamiltonian, defined by 

H T = H(Q)+ J £Q 2n+t Q l (Q) , (6) 



£ 



where Qi = Q n+ i = Pi, Q2n+e — v (i — 1,2, ■ • -,n). Note that H and f2 
only depend on Qi, ■ ■ ■, Q 2n . The non- vanishing a a 's are given by a, = Q n +i- 
The Euler-Lagrange equations of motion read 

= , (7) 

where is the Euler derivative 

= ®L - * ( -?£-) . (8) 
a dQ a dt \dQj 1 ) 

They follow from the variational principle 

5 f dt L(Q, Q) = £ / dt E®6Q a = , (9) 

where we have dropped a total derivative term in the last integral, and the 
5QaS are arbitrary variations. Local symmetries of the theory correspond to 



1 In the following we shall apply standard Lagrangean techniques to this first order 
Lagrangean. They are not to be confused with the symplectic approach discussed in [11]. 
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particular variations SQ a for which (9) vanishes without making use of the 
equations of motion. 

The merit of having recast the variational problem 5S T = in the above 
form, is that we can make use of well established Lagrangean methods [10] 
to extract possible local symmetries of the Hamilton equations of motion. In 
order to elucidate the basic ideas we shall consider in the following the case 
where we are dealing with a purely first class theory with only one primary 
constraint fl±. 

Consider the Euler derivative (8). From (5) one finds that it has the form 

e ( ^ = T,fS(Q)Q, + k^(Q) , (io) 

where (a, (3 = 1, • • •, 2n + 1) 

F^{Q) = d a ap-d p a a , (11) 

and 

H T = H(Q) + Q 2n+1 n 1 (Q) . (13) 

The only non-vanishing components of the (2n + 1) x (2n + 1) matrix 

are given by F^^ = —F^ +i — 1 (i — 1, • • •, n). Thus the matrix F^p has 
the form 

/ -1 \ 



I 

V o o o J 



(14) 



where 1 is an n x n unit matrix, is the n-component zero column vector, 
and is corresponding row vector. We now look for left-zero modes of F^ . 
The only zero mode is given by 

^ 0) = (0,0,1), (15) 

where is the n-component null vector. Hence 

^o) . £(o) = ^o) . £(o) = _ Qi (16) 
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On the zeroth level of the Lagrangean algorithm we therefore just recover the 
primary constraint. Per construction, it vanishes on the space of solutions to 
the Euler Lagrange equations of motion (7) (i.e., on shell). This constraint 
is only a function of the "coordinates" Q 1 , • • •, Q 2n - Hence its time derivative 
is given by 

ft ' (17) 

We next look for further constraints hidden in the equations of motion. To 
this effect we adjoin (17) to the equations of motion and consider the vector 

= ( -? ) = ( ) (18) 

which by construction vanishes again on-shell. The 2n + 2 components of 
can again be written in the form 

E S = E F a^iQ)Qp + K £(Q) 5 «i = l,2,---,2n + 2, (19) 

p 

where is the rectangular (2n + 2) x (2n + 1) matrix 

F« = ( ) (20) 

with xi the row-vector (-Vfti,0), and i?^) = (i?(°),0) = (-V// T ,0), where 

V = (a 1 ,,9 2 ,---,a 2n ) . (21) 

The zero modes of include the previous one, augmented by a zero (which 
will merely reproduce the zeroth- level constraint), and a new zero mode given 
by 

= (-a n+ ini,..., -^i,^!,...,^!, 0,1) . (22) 

One then readily verifies that 

vW . E w = ^ ■ K w = {H T , fix} , (23) 

where the Poisson bracket is taken with respect to q and p (In the following 
Poisson brackets are always understood to be taken with respect to q and p). 
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Since in the present case we only have one primary constraint, this Poisson 
bracket ist just {H, Oi}. It is only a function of the q^s and piS, i.e., of 
Qi, 4 ■ ■;Q2n- Let us assume, e.g., that it does not vanish on the primary 
surface. Then 

-0 2 = {F,0 1 } = ^-EW (24) 

is a secondary constraint, vanishing (by construction) on shell. We now 
proceed with the Lagrangean algorithm by adjoining the time derivative of 
^2 to the equations of motion. We are therefore led to consider the vector 



(25) 





( m \ 


i 


^(0) \ 


E^ = 




= A 

dt 


^(0) . #(0)) 




v ) 


\ dt 





E^ can again be written in the form 



(26) 



where a 2 = 1, 2, • • •, 2n + 3, and 



\ 




(27) 



with x.2 the row- vector (— Vf22,0). The non- vanishing components of the 
In + 3 component vector are again those of K^. F (2) possesses, besides 
the previous eigenvectors, augmented by an apropriate number of zeros, a 
new eigenvector having again the generic structure 



#( 2 ) = (-Q n+1 n 2 , ■ ■ -d 2n n 2 , <9i^ 2 , • • •, d n n 2 , o, o, i) 

Because of this generic structure, it follows that 

^■E^ = ^-K^ = {H T ,Q 2 } , 



or 



v 



< 2 > -E® = {H,Q 2 } + Q 2n+1 {n 1 ,Q 2 } 



(28) 

(29) 
(30) 



2 In order to demonstrate the general method we have to consider the case where at 
least one further constraint is generated. 
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Note that now the RHS depends on the Lagrange multiplier v = Q2n+i, if 
the Poisson bracket of the constraints (always taken with respect to q and p) 
does not vanish strongly. If we were to adjoin the time derivative of ■ E^ 
to the equations of motion, we would be led to an E^ 3 \ and corresponding 
rectangular matrix which possesses no further zero modes. A further 
constraint can however be generated if {Qi,Q 2 } vanishes on the subspace 
defined by fl 1 = Q 2 = 0, i.e., if 

{Q 1 ,Q 2 } = j2Ci 2 n e . (31) 

1=1 

In this case we can construct a further function of only Qi,-~, Q2n (i-e., which 
does not depend on the Lagrange multiplier v) , which vanishes on shell. This 
function is given by Q 

-n 3 = {H,n 2 } = ^ ■ - Q 2n+1 {n u n 2 } 

= ^■m+Q 2n+l j^C{ 2 {^-E^)), (32) 

t=\ 

where we have made use of (30), (16) and (24). If fi 3 is not a linear combina- 
tion of and ^2 then we construct E^ by adjoining the time derivative of 
fi 3 to the equations of motion. Proceeding in this way we (possibly) generate 
successively new constraints given by 

- n e+1 = 00 ■ m - Q 2n+l £ c k u n k . (33) 

k=l 

where the C\ t are defined by 

i 

{n 1 ,n e } = ^2c1 e (q,p)^k . (34) 
fc=i 

Note that the iterative process allows us to express the RHS of (33) in terms 
of scalar products ■ E^ with n < t. The algorithm will come to a halt 

3 Note that the constraints are generated precisely in the same systematic way as given 
by the Dirac algorithm. There exists therefore a complete correspondence between the 
Lagrangean and Dirac algorithm. 
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once no new constraint is generated. If v^ M ^ is the last zero mode generated 
by the iterative procedure, then, for a first class system, this will happen 
when v^ M "> ■ £?( M ) is a linear combination of the constraints fii, • • •, Vl M . We 
are thus left with a so-called "gauge identity", 

M 

0W.£W-£i^ = O . (35) 
e=i 

The Q-dependent coefficients K M can be related to the Poisson brackets 
{Qi,Qm} and {H,Qm} as follows. Because of the generic structure of the 
eigenvectors 

ti® = (-d n+1 n e , ■ ■ -a 2n n e , d^e, d n n e , o, • • •, o, i) ; i > 1 (36) 

we have that 

#M . EPO = {H T ,n M } , (37) 

or 

^M) . £(M) — ^jj Q M y + , ft M } . (38) 

Furthermore, the constraints generated by the Lagrangean algorithm are 
identical with those generated successively by the Dirac algorithm. If these 
constraints are first class, then, in particular, 

M 

{n 1 ,n M } = Y,ci M n e , (39) 

and (35) and (38) imply that 

M 

{H,n M } = J2h e M ^e ■ (40) 
1=1 

Hence 

K e M = h e M + vC[ M . (41) 

Note that in general h e M and C[ M depend on Q 1 , ■ • •, Q 2n - The constraints 
appearing in (35) can all be expressed in terms of ■ with t — 1, • • 
•,M — 1. It is therefore evident that the gauge identity (35) can be written 
in the form 

M-\ 

^M) . £(M) + £ p ( M )(Q)(^) . m) = . (42) 
e=o 
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Because of the iterative way in which the are computed, and the generic 
structure of the eigenvectors (36), this expression can be reduced step by 
step to the form 

2n+l M 

EEf(«) = 0, (43) 
a=l e=o 

where 

^ = *™{«*---^)- w 

The involve linear combinations of time derivatives up to (M — £)'th 
order. Multiplying the gauge identity (43) by an arbitrary function of time 
a(t), integrating the expression over time, and making a sufficient number 
of partial integrations (dropping total derivative terms, assuming that a(t) 
and time derivatives thereof up to M'th order vanish at the boundaries of 
integration), one is led to an identity of the form 



/ 



M /2n+l \ 

^EnE W (45) 

£=0 V a=l / 



The e £ 's depend on Q and time derivatives of Q, as well as on the arbi- 
trary parameter ct(t) and time derivatives thereof. For infinitessimal a(t) 
expression (45) has the form 

]T / dt Ef>8Q a = , (46) 

where 

M 

SQ a = -Y,^. (47) 
e=o 

The minus sign has been included for convenience. The integral (46) is just 
the variation of the action induced by the transformation Q a — > Q a + 5Q a 
(see eq.(9)). But because of the generic structure of the eigenvectors (36), 
we have that (recall that Qi = qi, Q n +i = Pi, Qin+\ — v ', i — 1, 2, • • •, n) 

M M 

i = i oqi f= . 
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Sv = -e° 



(48) 



These expressions have precisely the form following from the Dirac conjec- 
ture. In fact, the integral (45) can also be cast into a Hamiltonian version. 
To this effect we notice that it follows from (10), (12) and (14), and the 
structure of the eigenvectors (36) that 



2n+l 



a=l 



IF 



+ {H T ,n e } ; £> 1 . 



(49) 



Hence (45) can also be written in the form 



dt $V (p e + {H T , Q e }) + Ml, 

J U=i 



o, 



(50) 



where we have also made use of the identification e° = — Sv. This is pre- 
cisely the form of the variation of the action (4), with the variations Sqi and 
Spi given by (48). Indeed, for the case of only one primary constraint, the 
variation (4), written in these variables, has the form 



dt 



. dpi dpi J 



. dH XlA ' 



= 0, 
(51) 

which, for Sqi, Spi, and Sv given by (48), then becomes the following identity: 

J dt (^tot - t £ {H, n e } - v(?{n 1} n e } - Svn^j = o . (52) 

This is nothing but Eq. (50), after making a partial integration. If we are 
dealing with a first class system, then 



{if, Q e } = ^2 hfak , 

k 



Thus (52) becomes 



J dt e e {h{, + vCip^j ft e - SvSli) = . 



(53) 
(54) 

(55) 
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Hence the e e, s and 5v, determined by the Lagrangean algorithm, should be 
solutions to 

j a 

— -J(ht + vCb) = 0l (a = 2, 3), (56) 

Sv = — - e e (h\ + vCl e ) . (57) 

These equations can be solved iteratively [7,12]. The solution involves one 
arbitrary function of time (corresponding to one primary constraint). 

We now verify explicitely, for the case of an arbitrary first class system 
with one primary and two secondary constraints, that the parameters e e in 
(45) determined by the Lagrangean algorithm, are indeed solutions to Eqs. 
(56). For such a system the gauge identity will be generated at the third 
level of the Lagrangean algorithm. Independent of the specific form of the 
Hamiltonian one finds, proceeding in the manner described above, that the 
gauge identity (42) reads 



tfC). £(3) + R i^o) . ft )) + Ki&V -E 



+ K 



^).E^+vJ2C( 2 (^-E^) 



(58) 



where is given by (41) with M — 3. Making use of the generic structure 
of the eigenvectors (36), and of the recursive construction of one then 
verifies after some algebra that the gauge identity can be recursively reduced 
to the following form, 

« + (tfJ + |)(.W) 

+ j f v(Cl 2 + C? 2 |)K«i 0) 4 0) ) = , (59) 

where a summation over a = l,2,---,7is understood. The time derivatives 
act all the way to the right. By multiplying this expression with an arbitrary 
function of time e(t), and integrating the expression over time, dropping a 
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total derivative, we find, upon making use of (49), that it can be written in 
the form (50), where 

e 3 = e(t), 

e 2 = Kle-e, 

e 1 = e-±(Kte)-vC 2 12 e + Kle + vKtC 2 12 e. 

Sv = e 1 + vC\ 2 e - K\e - vK\C\ 2 e . (60) 

One can now easily verify that (60) are indeed solutions to (56) and (57), 
i.e., to 

e a -]Te% a = 0; (a = 2,3) , (61) 

3 

e 1 - ]T e l K\ - Sv = . (62) 
e=i 

Setting e 3 = e(t), and making use of the fact that, because of the way the 
constraints have been generated by the Lagrangean algorithm, h\ = Cf± = 
Cf 2 = 0, and h\ = h 2 = — 1 in (56) and (57), one can solve the above 
equations recursively, starting with equation (61) for a = 3. The solutions 
are given, as expected, by (60). 

Concluding, we have shown that well established Lagrangean techniques 
can be used to uncover all the gauge symmetries of the Hamiltonian equations 
of motion for first class systems with one primary constraint, and that there 
is a direct correspondence between each level of the Lagrangean algorithm 
and the Dirac approach to local symmetries based on the Dirac conjecture. 
This analysis can be extended to systems with several primary constraints. 
A general discussion of gauge symmetries for arbitrary systems based on the 
approach presented in this paper, is deferred to a future publication. 
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